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Abstract 


Yubin, Z., Zhengying, W., Lei, Z., Qinyin, L., £ Jun, D. 
(March-April, 2017). Improved online sequential extreme 
learning machine for simulation of daily reference 
evapotranspiration. Water Technology «and Sciences (in 
Spanish), 8(2), 127-140. 


The traditional extreme learning machine has significant 
disadvantages, including slow training, difficulty in 
selecting parameters, and difficulty in setting the singularity 
and the data sample. A prediction model of an improved 
Online Sequential Extreme Learning Machine (IOS-ELM) 
of daily reference crop evapotranspiration is therefore 
examined in this paper. The different manipulation of the 
inverse of the matrix is made according to the optimal 
solution and using a regularization factor at the same time in 
the model. The flexibility of the IOS-ELM in ETO modeling 
was assessed using the original meteorological data (Tmax, 
Tm, Tmin, n, Uh, RHm, Y, Z) of the years 1971-2014 in 
Yulin, Ankang, Hanzhong, and Xi'an of Shaanxi, China. 
Those eight parameters were used as the input, while the 
reference evapotranspiration values were the output. In 
addition, the ELM, LSSVM, Hargreaves, Priestley-Taylor, Mc 
Cloud and IOS-ELM models were tested against the FAO- 
56 PM model by the performance criteria. The experimental 
results demonstrate that the performance of IOS-ELM was 
better than the ELM and LSSVM and significantly better 
than the other empirical models. Furthermore, when the 
total ETO estimation of the models was compared by the 
relative error, the results of the intelligent algorithms were 
better than empirical models at rates lower than 5%, but 
the gross ETO empirical models mainly had 12% to 64.60% 
relative error. This research could provide a reference to 
accurate ETO estimation by meteorological data and give 
accurate predictions of crop water requirements, resulting in 
intelligent irrigation decisions in Shaanxi. 


Keywords: Daily reference evapotranspiration, extreme 
learning machine, online learning, matrix singularity. 


Resumen 


Yubin, Z., Zhengying, W., Lei, Z., Qinyin, L.,6 Jun, D. (marzo- 
abril, 2017). Máquina de aprendizaje extremo secuencial en línea 
mejorada para la simulación de la evapotranspiración de referencia 
diaria. Tecnología y Ciencias del Agua, 8(2), 127-140. 


La máquina de aprendizaje extremo tradicional tiene desventajas 
significativas, tales como entrenamiento lento, dificultad en la 
selección de parámetros y dificultad en establecer la singularidad y 
la muestra de datos. Por lo tanto, en el presente artículo se examina 
un modelo de predicción de una máquina de aprendizaje extremo 
secuencial en línea mejorada (1IOS-ELM) de la evapotranspiración de 
referencia diaria de cultivos. La diferente manipulación de la inversa 
de la matriz se hace de acuerdo con la solución óptima y utilizando un 
factor de regularización al mismo tiempo en el modelo. La flexibilidad 
de la 1OS-ELM en la modelación de la ETO se evaluó empleando los 
datos meteorológicos originales (Tmax, Tm, Tmin, n, Uh, RHm, q, 
Z) de los años 1971-2014 en Yulin, Ankang, Hanzhong, y Xi'an en 
Shaanxi, China. Estos ocho parámetros se usaron como insumos o 
datos de entrada, mientras que los valores de la evapotranspiración 
de referencia fueron los datos de salida o el producto. Asimismo, se 
probaron los modelos ELM, LSSVM, Hargreaves, Priestley-Taylor, 
Mec Cloud y 1OS-ELM contra el modelo FAO-56 PM mediante los 
criterios de desempeño. Los resultados experimentales demuestran 
que el desempeño de IOS-ELM fue mejor que le de ELM y LSSVM 
y significativamente mejor que los demás modelos empíricos. Más 
aún, al comparar la estimación total de ETO de los modelos mediante 
el error relativo, los resultados de los algoritmos inteligentes fueron 
mejores que los modelos empíricos a índices inferiores a 5%, pero 
los modelos empíricos de ETO bruta tuvieron un error relativo de 
12 a 64.60%. Esta investigación podría proporcionar una referencia 
para la estimación precisa de ETO mediante datos meteorológicos y 
proporcionar predicciones precisas de los requerimientos de agua de 
los cultivos, lo cual resultaría en decisiones de riego inteligentes en 
Shaanxi. 


Palabras clave: evapotranspiración de referencia diaria, máquina 
de aprendizaje extremo, aprendizaje en línea, singularidad de la 
matriz. 
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Introduction 


The calculation of reference crop evapo- 
transpiration is a key to intelligent irrigation 
systems. Therefore, accurate estimation of ETO 
becomes important in irrigation schedules for 
planning and optimizing the agriculture area. 
Numerous methods have been put forward to 
estimate ETO. Of these, the Penman-Monteith 
56 (PM) model has given the best results; it 
was officially recommended by the Food and 
Agriculture Organization of the United Nations 
(FAO) in 1998 (Allen et al., 1998). FAO selected 
the PM model as the standard equation for ETO 
estimation because it can provide the most ac- 
curate results in the world. ETO is also a key 
parameter in the design of intelligent irrigation 
for field crops. Engineers need to know the ir- 
rigation water consumption requirements for 
each crop so that they can calculate or estimate 
the remaining components of the water balance. 
Also, agriculturists need to obtain the specific 
water requirements of a crop so that they can 
generate a satisfactory yield. Itis also necessary 
to know whether these specific requirements are 
being met with ordinary irrigation. 

As described by Kisi (2008) and (Yubin et al., 
2014), about 50 measures have been proposed 
for estimating evapotranspiration, which can be 
sorted into four types: radiation, temperature, 
synthetic and evaporating dish. FAO assumed 
the ET definition given by Smith, Allen and 
Pereira (1997), and adopted the FAO-56 PM as 
the sole equation for estimation of ETO. 

However, the PM model requires a lot of 
meteorological data as input, and the calculation 
process involves complex and nonlinear regres- 
sion among these factors. So, a simpler and 
more accurate simulation model needs to seek 
ETO in the case of lack of meteorological data. 
Thus, artificial intelligence algorithms (e.g., neu- 
ral networks) for reference evapotranspiration 
(ETO0) modeling have been given more attention 
in recent decades. Feng and Cui (2015) found 
that an ELM model gave better results than 
empirical models in the area of central Sichuan. 
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Kisi (2007) estimated daily ETO using the ANN 
method and compared their calculation results 
with the other models. Ozgur Kisiet (2013) pro- 
posed a reference evapotranspiration model by 
LSSVM. Kisi (2011a) considered daily ETO using 
wavelet regression model and compared this 
model to other empirical models. Kisi (2011b) 
modeled ETO using evolutionary feed-forward 
neural networks. Marti, Gonzalez-Altozano and 
Gasque (2011) used ANN to estimate daily ETO 
without local climatic data. Kumar, Raghuwan- 
shi and Singh (2011) researched the application 
of ANN in estimating evapotranspiration 
modeling. Shiri et al. (2012) established an ETO 
simulation model using GEP (gene expression 
programming) for Spanish Basque, and found 
that the GEP model performed better than the 
ANFIS, Hargreaves and Priestley-Taylor mod- 
els. Wang, Traore and Kerh (2008), and Traore, 
Wang and Kerh (2010) estimated daily ETO 
using BP-ANN. However, BP-ANN has major 
disadvantages, such as its slow speed of training 
and difficulty in selecting parameters. In recent 
years, new intelligent algorithms have appeared 
in the industrial field, such as extreme learning 
machines (ELM) and support vector machine, 
among others. 

In the present study, ELM is proposed as 
an alternative to other models for predictive 
control. It can randomly choose hidden nodes 
and analytically determine the output weights 
of SLENs. However, ELM cannot confirm the 
singularity of the output matrix of the hidden 
layer, and it also cannot make fine tuning ac- 
cording to the characteristics of the data set, 
which will affect its efficiency and stability. 

The main purpose of this paper is to opti- 
mize the ELM approach in the modeling of daily 
ETO using the original meteorological data. All 
previous studies have indicated that intelligent 
models can input the factors of FAO-56 PM as 
they estimate ETO. In fact, these factors will be 
another complex computing project by meteoro- 
logical raw data, to avoid creating more severe 
error during multistage formula calculation. 
Also, the manipulation of the inverse of the 
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matrix is adjusted with reference to the optimal 
solution and the regularization factor at the 
same time, which is motivated by the online 
learning method. In summary, an improved 
online sequential extreme learning machine 
(IOS-ELM) is designed, and the new algorithm 
can produce good generalization performance 
in a model of daily ETÚ in an irrigation system. 


Materials and case study 


The study was conducted in Yulin (38.27? 
N, 109.78? E), Ankang (32.72% N, 109.03? E), 
Hanzhony (33.07%N, 107.7*E) and Xi/an (34.39 N, 
108.9? E) in Shaanxi province in China, shown 
in Fig. 1. The area has a hot and dry climate for 
the greater part of the year. 

Daily meteorological data used for this 
study was from the years 1971-2014. The fol- 
lowing observed eight meteorological variables 
with daily temporal resolution were used: 
wind speed at 10m above the ground (Uh), 
mean temperature (7), mean relative humidity 
(RH), minimum temperature (Tmin), maximum 
temperature (Tmax), actual Sunshine duration 
(n), latitude ($) and elevation (Z), which were 
downloaded from China meteorological data 
sharing service system (http: / /cdc.nmic.cn/ 
home.do). Data from the first 29 years (1971- 
1999) was used to train the models. Data from 
the next ten years (2000-2009) was used for the 
test. The data from the remaining years was 
used for validation. It must be noted, however, 
that missing data was replaced by the average 
of the data from the day before and the day 
before. The regional climate characteristics are 
given in table 1. 


Figure 1. The location of the cities. 


Methodology 


Calculation models of reference crop 
evapotranspiration 


The study focused on the comparison of the pro- 
posed IOS-ELM model with the ELM, LSSVM, 


Table 1. Means of main variables. 


Uh T RH T min T max n (0) YA 
Yulin 2.14 9.2 53.5 3.1 16.2 7.32 0.67 1157 
An kang 1.35 15.9 74.1 12.1 21.4 4.58 0.57 290.8 
Han zhong 1.15 15.1 78.5 11.5 19.8 4.03 0.58 509.5 
Xian 1.60 14.6 64.3 9.7 19.6 4,54 0.60 397.5 
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Hargreaves, Mc-Cloud, and Priestley-Taylor 
models. First, the ETO values of four cities were 
calculated using the FAO-56 PM. Then, the stan- 
dard formula evapotranspiration calculation for 
all empirical models is shown. 


(1) FAO-56 Penman-Monteith: 


0.418A(R, -G)+ a (e, -e,) 


ETO- 
A+y(1+0.34u,) 


(1) 


The original meteorological data of T, 


T, T.,, 1, Uh, RH,,, $ and Z were used in the 


model: 


(2) Hargreaves: 


0.0023xR, (7+17.8)(T,,,. -T,, 
A 


ETO = 


T 


max 


T,T 


min 


n and q were used in the model: 
(3) Mc-Cloud: 
ETO =0.254 x 1.079" (3) 
Only T was referred to in the model: 


(4) Priestley-Taylor 


ero=1262_.[(Rr-6) (4) 
A+y A 


T ,T,T 


min 


e n and q were used. 

However these variables are obtained 
directly or indirectly from the meteorological 
raw data (Uh, T, RH, Tmin, Tmax, n, q and Z). 
Furthermore, the calculation formula for them 
did not have a precise formula by estimation or 
experience. 

Therefore, the inputs Uh, T, RH, Tmin, Tmax, 
n, $ and Z, the ETO output were calculated 
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by the FAO-56 PM method and used for the 
calibration of the IOS-ELM models. The mean 
absolute error (MAE), the root mean square 
error (RMSE), effectiveness index of the model 
(EF) and self-correlation coefficient (R?) statistics 
were used for the assessment criteria of the 
models in this study. EF model efficiency mainly 
depends on the Nash coefficient EF values; as 
the values approach one, the efficiency of the 
model increases. The study adopted the calcula- 
tion model of the validity index for EF by Nash 
and Sutcliffe. 


Extreme learning machine (ELM) 


For N random distinct samples (x, t) where x, 
rl ER ds lo ER” and 


for the standard SLENs ( Ñ hidden nodes), the 
activation function g(x) is expressed as: 


j=1,...,N (5) 


where vw, = [w,,, 0,,,..., 1, ]'is the weight vector 
connecting the ith hidden node and the input 
nodes, P,=[B,,, B;,-.., B, ]'is the weight vector 
connecting the ith hidden node and the output 
nodes and b, is the threshold of the ith hidden 
node. 


pg (w,-2,+b)=+, j=1,...,N (6) 
i=1 


The above N equations can be written com- 
pactly as 


Hf=T (7) 


where 


T T 
1 t, 

p=IM| ,T=IM (9) 
Bs Ñxm y Nxm 
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where H is called the hidden layer output ma- 
trix in the neural network, and the ith column of 
H is the ith hidden node output toward inputs 
Mid 

The aim is to solve the above issues and 
put forward an extreme learning machine for 
SLENSs. 

A training set was provided as: 


N 
N = (8 4), ER, e r") , the activation 
function g(x), and hidden node umbel Ñ. 
Step 1: Randomly allocate input weight v, 
and bias b, ¡=1,2,..., Ñ. 
Step 2: Calculate the hidden layer output 
matrix H. 


Step 3: Calculate the output weight f. 
$p=HT (10) 


where T=[f,, ..., t, 1*, H' is a generalized inverse 
of MP. 


Online sequential ELM (OS-ELM) 


ELM is a relatively effective and simple al- 
gorithm that is also able to learn quickly and 
generalize well. However, meteorological data 
are difficult to collect and the data set is large, 
which may cause a decline in the performance 
of the ETO model. Thus, the online sequential 
extreme learning machine (OS-ELM) by Liang 
(2006) was referenced in the previous research. 

The output weight matrix f = H"Tis a least- 
squares solution of (7). Meanwhile, the matter 
where rank(H) = Ñ the number of hidden nodes 
(Ao, Xiao, €: Mao, 2009) is considered. So, H*of 
(10) is given as: 


-1 
(Ha) H (11) 
Tf H7H tends to become fantastic, it can also 
be made nonsingular by increasing the number 


of data or choosing a smaller network size. 
Substituting (11) into (10) gives: 


p= (474) HT (2) 


Equation (12) is called the least-squares 
solution to Hf = T. Sequential implementation 
of the least-squares solution of (12) gives the OS- 
ELM. However, the OS-ELM may have some 
deficiencies, especially the fact that solving the 
generalized inverse matrix MP of H may cost a 
huge amount of time in the training process. The 
general method of singular value decomposition 
is used to solve matrix H, but its computational 
complexity is O(4NNÑ? + 8N*) (Brown, 2009). 


Improved algorithm of OS-ELM (10S- 
ELM) 


This paper proposes an improved OS - ELM 
called IOS-ELM. This new model was devel- 
oped by modifying and improving the singu- 
larity of the matrix. First, Equation HB = T 
will be replaced by H'Hf = H"T, which has at 
least one optimization solution. This reduces 
the computational complexity of solving the 
inverse, which results in a reduction of the train- 
ing time. Second, the regular factor 1/2 is joined 
when calculating the output weights. Last, the 
subsequent online learning stage is added. In 
theory, this algorithm can provide good gen- 
eralization performance at an extremely fast 
learning speed. 

Step (1): Allocate random input weights vw, 
and bias b, initialize network and calculate the 
initially hidden layer output matrix H,. 

Step (2): Set r = rank(H), if r = N,, then cal- 
culate the initial weight matrix $ =P HT, 1£ 
r = N then calculate the initial weight 
p?= HIP, T.. 


AE -1 
Where P, > EH) ¿E (EH) . 


If r « N, and r « Ñ, to solve the two optimiza- 
tion models: 


min|M-B[ and min [Bp- Cc 


BEg BER! 


Then, the optimization solution B* and f' can 
be obtained. 
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Where c= HIT, M= HH, 3*=(g ERYo No), 
gis a positive definite symmetric matrix. 

Step (3) Set K= 0; then, present the (K + 1)th 
chunk of new observations: 


where N,, is the number of observations in the 
(K + Dth chunk. 

Step (4) Calculate the partially hidden layer 
output matrix H, , for the (K + 1)th chunk of 


data N,.., as shown in (17): 


K+1/ 


3 (%,:x,+b,) - g(i0,-x,+b,] 
Ho = : Es : (13) 
g(10,'xy,, +b,) .. g(1,:x,+b,] 


N q 


Step (5) According to step (1), calculate the 
output weight P**1, 
Step (6) Set k =k + 1. Go to Step (3). 


Application and results 
IOS -ELM model under lack of data 


The original eight meteorological parameters 
chose and combined with a different pattern 
in this section, which was taken as input 
values. Meanwhile, the calculation of FAO 
56 Penman-Montieth was put as the output 
value. By this method, the IOS-ELM model 
is established. However, it needs to further 
consider the effectiveness of the combination 
pattern among the eight meteorological data. 


Therefore, the correlation between ETO0 and the 
data was analyzed. In this way, ISO-ELM can 
choose reasonable meteorological parameters to 
complete the forecast even if there is a lack of 
meteorological data. This is shown in table 4. 

It can be clearly seen from table 2 that the 
ETO outperformed all eight meteorological pa- 
rameters in terms of correlation. Although the 
data set is not similar for different cities, the cor- 
relation behaved in the same way. Tis closely 
correlated with evapotranspiration for each city, 
followed by the average temperature, minimum 
temperature, the actual sunshine time and wind 
speed. The influence of the latitude and altitude 
were so small that they were negligible. Finally, 
the humidity is negative. 

The simulation accuracy of the TOS-ELM 
model was discussed by referring to table 2 
under lack of meteorological data. It should 
be noted that the latitude and altitude were 
eliminated because they had virtually no effect 
on the results. Taking Yulin city as an example, 
the first 10-year (1971-1999) span of data was 
used to train the models. Then, using different 
combinations, the error was analyzed, as well 
as the correlation coefficient and effectiveness 
of the prediction. This is shown in table 3. 

The ISO-ELM model was applied by com- 
paring the different parameters shown in table 
3. It is immediately noticeable that the predic- 
tion results were the same for eight-parameter 
and six-parameter inputs. That is because the 
latitude and altitude almost have no effect on 
the prediction for the same station. Secondly, 
the temperature had the largest influence on 
the prediction, particularly the maximum tem- 
perature. As long as the temperature is one of 


Table 2. Correlation of data and ETO. 


Uh mí RH,, Le 1d " $ z 
Yulin 0.20 0.88 -0.36 0.79 0.92 0.50 -2e-15 
Ankang 0.15 0.85 -0.34 0.71 0.93 0.70 -7e-16 ie 
Hanzhong 0.21 0.85 -0.48 0.73 0.93 0.69 de-15 
Xian 0.30 0.80 -0.43 0.72 0.85 0.66 Le-16 
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Table 3. Influence of different meteorological data combinations on ETO. 


Model inputs RMSE R? EF% 

All 0.4132 0.9696 95.72 
O Ra 0.4132 0.9696 95.72 
Duo E Tete 0.7865 0.9625 92.3 
O dl 0.6737 0.9619 94,3 
T,, T,T,., Uh, RH 0.7868 0.9620 92.2 
Tias 1, 1, Uh, RH 0.9021 0.9496 89.8 
To. Ty 1 Uh, RH 0.9135 0.9490 89.6 
T, T,..., 1, Uh, RH 0.9666 0.9423 88.3 
os Y Ty te 0.8270 0.9591 91.4 
Tay IE blo 0.7338 0.9673 93.2 
E, a O 0.6368 0.9651 94.9 
T,..., 1, Uh, RH 1.2321 0.9060 81.1 
T,n, Uh, RH 1.0048 87.4 
T,..., 1, Uh, RH 0.8428 91.1 
ON 0.8505 90.9 
n, Uh, RH 2.2612 36.3 
os RE 0.6480 94.7 
IA 0.8580 90.8 
A 0.7942 92.1 
T,., Uh, RH 0.7436 93.1 
Tn, RH 0.6298 95.06 
TE 107 WTO 0.8786 90.3 
ES 0.9464 88.8 
a 0.9119 89.64 
Tas RH 0.6258 95.12 
EA 0.9012 0.9483 89.89 
.. y Uh 0.9999 93.59 87.55 
n, RH 2.1127 0.6694 44.42 
n, Uh 2.1611 0.6543 41.85 
Uh, RH 2.6964 0.3395 9.47 


the parameters, the model is accurate. Thirdly, models including Hargreaves, Mc-Cloud and 


when only two temperatures were used as the 
inputs, it still performed better than RMSE, EF 
and R?statistics. So, properly reducing some 
variables and adopting reasonable combina- 
tions of variables can improve the accuracy of 
prediction. 


Comparison with other calculation 
formulas 


The IOS-ELM model was compared with the 
ELM and LSSVM, as well as conventional 


Priestley-Taylor methods in respect of RMSE 
and MAE statistics in different cities in tables 
4-5. There are six parameters as input variables 
in the model. 

Tables 4-5 show that IOS-ELM outperformed 
all other models by all performance criteria. 
Compared with the intelligent and empirical 
models, the ISO-ELM performed the best value 
of RMSE<0.46 and MAE<0.41, and the ELM 
and LSSVM models performed better than the 
others. A few differences appeared among the 
Mc Cloud and Priestley-Taylor models. It was 
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Table 4. RMSE of the models in the test period. 


A RMSE (mm/day) 
Yulin A kang Hanzhong Xian 
IOS-ELM 041 0.45 0.45 0.41 
ELM 0.86 0.72 0.88 0.78 
LSSVM 0.96 1.38 1.30 1.05 
Hargreaves 2.18 1.38 127 0.79 
Mc Cloud 3.53 2.20 1.98 1.70 
Priestley-Taylor 2.43 1.73 1.59 0.71 
Table 5. MAE of the models in the test period. 
AN MAE (mm/day) 
Yulin Ankang Hanzhong Xi'an 
IOS-ELM 0.40 0.35 0.33 0.31 
ELM 0.52 0.68 0.62 0.55 
LSSVM 0.77 1.20 1.13 0.93 
Hargreaves 1.97 1.25 1.19 0.58 
Mc Cloud 3.15 1.86 1.66 1.36 
Priestley-Taylor 2.13 1.56 1.40 0.51 


also discovered that the Hargreaves method 
provided better accuracy than other methods 
among the empirical models. 

Although IOS-ELM, ELM and LSSVM mod- 
els had better simulation effects, the running 
time is distinguishing, as shown in table 6. 

It is clear from table 6 that the TOS-ELM 
model runs faster than ELM and LSSVM in the 
process of calculating by at least 24.8%. 

In order to consider the portability and error 
causes of the IOS - ELM model, the estimates of 
each model for four cities are shown in figures 
2-5 in the form of scatter plots in the validation 
period. It is generally clear from the scatter 
plots that the six input ISO-ELM estimates are 
closer to the corresponding FAO-56 PM ETO 


values than other models. The fit line equations 
y = ax + b and R? values indicate that the ISO- 
ELM model performed with better accuracy. 
Meanwhile, the a and b coefficients of the six- 
input ISO-ELM model were closer to 1 and 0, 
respectively, with a higher R? value than those 
of the other models. 

For Yulin, ISO-ELM and ELM estimates were 
closer to the FAO-56 PM ETO0 values than those 
of the other models (R?> 0.96). A slight differ- 
ence exists between LSSVM, and Hargreaves 
was better than the surplus models. The Mc 
Cloud estimate had the least accuracy. It can be 
concluded that the ISO-ELM and ELM models 
are the best methods to use for daily ETO estima- 
tion in Yulin. 


Table 6. Running time of different models. 


Running time (S) 
Model 
Yulin Ankang Hangzhong Xi'an 
IOS-ELM 16.5 17.9 13.8 17.1 
ELM 28.8 23.8 21.2 25.8 
LSSVM 22.4 30.5 19.9 36.5 
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For Ankang, ISO-ELM and Hargreaves were Mc Cloud estimate was also the least accurate 3 
closer to the FAO-56 PM ETO values, a slight (R?= 0.6141). This leads to the conclusion that in 3 
difference exists between LSSVM, and ELM was this city, the ISO-ELM and empirical Hargreaves 3 
. . a 
better than the Mc Cloud model. In this city, the models were the best. E 
ES 
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Figure 3. The FAO-56 PM ETO and estimated ETO values of Ankang. 


For Hanzhong, ISO-ELM was closer to the 
FAO-56 PM ETO of R?= 0.9911, followed by the 
Hargreaves, ELM, and LSSVM models. The Mc 
Cloud and Priestley-Taylor estimates were the 
least accurate. 
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For Hanzhong, ISO-ELM was closer to the 
FAO-56 PM ETO with R?=0.9905, followed by the 
ELM, LSSVM, Priestley-Taylor, Hargreaves and 
Mc Cloud models, which had R? values of 0.9547, 
0.9488, 0.8804, 0.8335 and 0.5577, respectively. 
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Figure 4. The FAO-56 PM ETO and estimated ETO values of Hanzhong. 


The total ETO estimation of every model is 
compared in table 7 because of its importance 
in irrigation management. The ISO-ELM clearly 
performed better than the other models from 
the relative error, which was 4.76, 0.23, 0.02 and 


0.54%, respectively. In four of the cities, it gave 
the closest estimate to the total FAO-56 PM ETO 
during the validation period. 

For Yulin, the ELM and LSSVM had the 
same accuracy, which was the second best, and 
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Figure 5. The FAO-56 PM ETO and estimated ETO values of Xian. 


had 5274 and 5488 estimates lower than the 
10% error, respectively. They were followed 
by Hargreaves, Priestley-Taylor and, lastly, 
McCloud (1955), which had the highest error 
at 64.60%. For Ankang, the LSSVM was ranked 
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as the second best, followed by the ELM, Mc 
Cloud, Hargreaves and, finally, Priestley-Taylor. 
For Hangzhong, the LSSVM was ranked as the 
second best, followed by ELM, Hargreaves, Mc 
Cloud, and Priestley-Taylor, which had 28.49, 
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Table 7. Performance statistics of the models in the validation period. 


Total ET0(mm) Relative error (%) 
Models 

Yulin Ankang | Hanzhong Xi'an Yulin Ankang | Hanzhong Xian 
Observed 5 039 4 661 4 307 3134 = = = = 
ISO-ELM 5279 4 650 4.306 3117 4.76 0.23 0.02 0.54 
ELM 5274 4 636 4413 3 207 4.66 0.54 2.46 2.33 
LSSVM 5 488 4.639 4. 408 3226 8.91 0.47 2.35 2.94 
Hargreaves 3 003 3 362 3 080 3276 40.40 27.87 28.49 4.53 
Mc Cloud 1784 3672 20 2.436 64.60 PAL 22 30.37 22.27 
Priestley-Taylor 2 844 3051 2864 2.727 43.56 34.54 33.50 12.99 


30.37 and 33.50% error, respectively. For Xian, 
the ELM, LSSVM, and Hargreaves were ranked 
the second best, followed Priestley-Taylor and, 
finally, Mc Cloud. 

In short, for the different cities, the ISO-ELM 
performed better than the other models, and the 
other models had different degrees to adapt to 
the application. 


Conclusion 


The improved sequential extreme learning ma- 
chine (IOS-ELM) is designed and applied for 
simulation of daily reference evapotranspiration 
through different manipulation of the inverse of 
the matrix and using the regularization factor 
and online learning method at the same time 
Experimental results demonstrated that the 
IOS-ELM can learn faster and achieve better 
performance than traditional ELM. 

First, the IOS-ELM model effectively over- 
comes the defects of traditional ELM, such as 
slow training speed, difficult parameter deci- 
sions, difficulty in setting the singularity and 
effect of data samples. 

Second, the potential of the ISO-ELM 


pirical models were limited to specific regions 
and required modification. 

Fourth, in the different meteorological data 
combinations for ET'O estimation, as long as 
there was a temperature-related parameter 
calculation, the calculation accuracy of ETO was 
over 94%, and Tmax was especially effective. 
These accurate calculations can be a valuable 
reference for the development of intelligent ir- 
rigation in water decision-making systems. 
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Notation 
The following symbols are used in this paper: 


ETO = reference evapotranspiration (mm day”). 

A  =slope of the saturation vapor pressure 
function (kPaC”). 

Rn  = net radiation (MJ m? day”). 

G  = soil heat flux density (MJ m? day”). 


technique for the estimation of reference c = psychometric constant (kPa C”). 
evapotranspiration was investigated for four T  = mean air temperature (*C). 

areas in Shaanxi of China; particularly, eight U2 = average24-h wind speed at 2 m height 
meteorological data were used as inputs. (ms”). 


Third, it was demonstrated that intelligent 
algorithm models (IOS-ELM, ELM, and LSSVM) 
are widely applicable to different areas, but em- 


Rs  = solar radiation (MJ m? day”). 
es 


the saturation vapor pressure (kPa). 
ea  = the actual vapor pressure (kPa). 
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